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Abstract 



Energy eigenvalues and order parameters are calculated by exact diagonal- 
ization for the transverse Ising model on square lattices of up to 6x6 sites. 
Finite-size scaling is used to estimate the critical parameters of the model, con- 
firming universality with the three-dimensional classical Ising model. Critical 
amplitudes are also estimated for both the energy gap and the ground-state 
energy. 
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I. INTRODUCTION 



Recent advances in computer technology have allowed the exact diagonalization of Ising- 
typc quantum spin systems up to 36 sites in size. Schulz, Ziman and Poilblanc (1996), 
for example, studied the J1-J2 XXZ Heisenberg spin model on square lattices up to 6x6 
sites. Our aim in this paper is to carry out an exact diagonalization study of the transverse 
Ising model on the square lattice, in order to estimate its critical parameters and study its 
finite-size scaling behaviour. 

The transverse Ising model in (2+l)D is well-known to be the quantum Hamiltonian 
corresponding to the classical 3D Ising model (Suzuki 1976, Pradkin and Susskind 1978), 
and exhibits a quantum phase transition in the same universality class as the classical 3D 
Ising thermal transition. It was first studied by series expansion methods by Pfcuty and 
Elliott (1971), and there have been several further series expansion calculations since then, 
both 'low-temperature' (Marland 1981, Yanase et al 1976, Oitmaa et al 1991) and 'high- 
temperature' (Hamer and Irving 1984, Hamer and Guttmann 1989, He et al 1990). Exact 
finite-lattice calculations have also been carried out prcvioTisly (Roomany and Wyld 1980, 
Hamer 1983, Hcnkcl 1984, 1987) for square lattices of up to 5x5 sites, and similar calculations 
have also been done for the triangular lattice (Hamer and Johnson 1986, Henkel 1990, Price 
et al 1993). Here we extend these calculations for the first time to the 6x6 lattice, and use 
finite-size scaling theory to obtain improved estimates of the critical point and critical index 

V. 

The finite-size scaling amplitudes at the critical point are also of interest. In (1 + 1)D, 
it is well-known that the theory of conformal invariance relates the scaling amplitudes to 
fundamental parameters of the underlying effective field theory at the critical point, such 
as the conformal anomaly and scaling indices. In higher dimensions, a similar scenario is 
known to hold at 'first-order' transitions, where a continuous symmetry is spontaneously 
broken, giving rise to Goldstone bosons (Hasenfratz and Leutwyler 1990): the finite-size 
scaling amplitudes are related to parameters of the Goldstone bosons such as the spin-wave 
stiffness and spin-wave velocity. Does something similar apply at second-order transitions 
in higher dimensions? Apart from a discussion by Cardy (1985), little has been done in 
this area. One peculiar result was obtained by Henkel (1986, 1987) and Weston (1990), and 
confirmed recently by Weigel and Janke (1999): the scaling amplitudes of the spin-spin and 
energy-energy correlation lengths on antiperiodic lattices have a universal ratio: 

4^ = — (1-1) 

where the Xi are the scaling indices in the respective sectors. This phenomenon appears to 
have no good theoretical explanation at the present time. 

Our exact diagonalization methods are outlined briefiy in Section 2, and the numerical 
results are presented in Section 3. The critical parameters so obtained are compared with 
other estimates in Section 4, and the critical amplitudes are discussed. 

II. METHOD 

The transverse Ising model on the square lattice has the Hamiltonian 
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i <ij> i 



(2.1) 



where the sum < ij > runs over nearest neighbour pairs on the lattice, and the a matrices 
are the usual Pauli spin operators acting on a 2-state spin- variable at each site. The coupling 
X is analogous to an inverse 'temperature', and h represents an external 'magnetic field'. We 
shall employ a representation in which the (T3{i) are diagonal. Periodic boundary conditions 
are assumed. 

The unperturbed ground-state of the model at x = has all spins 'up', i.e. (Jsii) = +1, 
all i. The interaction term will induce an admixture of states with 'flipped' spins. The 
Hilbert space of the model consists of two sectors, containing an odd and even number of 
flipped spins respectively. 

Exact diagonalizations have been carried out for LxL lattices, L = 1,..,6. The methods 
employed are fairly standard, for the most part, and will not be described in detail here. 
First, a list of allowed basis states in the given sector was prepared, using the 'sub-lattice 
coding' technique of Lin (1990). This efficient technique produces a sorted list of states, 
requiring only one integer word of storage per state. Since only the zero-momentum states 
are considered here, the states were 'symmetrized': that is, all copies of a given state under 
translations, refiections and rotations were represented by a single state. Thus for the 6x6 
lattice in the even sector, the total number of 'unsymmetrized' states is approximately 2^^, 
whereas under symmetrization this is reduced by a factor of approximately 288, down to 
119,539,680. 

Next, the Hamiltonian matrix elements are generated, by applying the interaction opera- 
tors of equation ( ^.1[ ) to each initial state, sjnnmetrizing the resulting final state, and looking 
it up in the master file. The elements were grouped into blocks, each of which acts between 
small sub-sets of the initial and final state vectors, to avoid 'thrashing' during the matrix 
multiplications. Within each sub-set, the initial and final addresses can each be fitted into 
a half-integer, so that the matrix elements occupied 35 Gbyte of storage over all. 

Finally, the lowest eigenvalue and eigenvector of the Hamiltonian were found in each 
sector, using the conjugate gradient method. Nightingale et al (1993) showed that the 
conjugate gradient method converges faster than the Lanczos method for large problems 
such as this. We found that the eigenvalue converged to an accuracy of 1 part in 10^° in 
20-25 iterations for the 6x6 lattice in the neighbourhood of the critical point. 

Having determined the quantities of interest for each finite lattice, it is then necessary 
to make an extrapolation to L ^ 00, to estimate the bulk behaviour of the system. In the 
vicinity of the critical point, the finite-lattice sequence will typically behave as 

/l = /oo + aiL--^ + a2L--^ + ■ ■ ■ (2.2) 

where the uji are non- integer exponents, in general (Barber 1983). The problem of extrapo- 
lating such a sequence has been discussed in several reviews (Smith and Ford 1982, Barber 
and Hamer 1982, Guttmann 1989). We have employed a number of different algorithms, 
including: 

i) the Neville table (Guttmann 1989), which is best suited to a simple polynomial se- 
quence, with integer exponents uf, 
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ii) the alternating VBS algorithm (van den Broeck and Schwartz 1979, Barber and Hamer 
1982), which can give good convergence for sequences of type (|2.2| ), but needs at least 
two iterations to work well; 

iii) the Lubkin algorithm (1952), which is more suitable for short sequences; 

iv) the Bulirsch-Stoer algorithm (1964), which has been applied in this context by Henkel 
and Patkos (1987), and Henkel and Schiitz (1988). This algorithm involves an explicit 
parameter uo which can be optimized to match the leading power-law correction. It 
has been claimed by Henkel and Schiitz that the algorithm is more robust and more 
accurate than the VBS algorithm, especially for short sequences. 



III. RESULTS 

A. Finite-lattice data 

The pseudo-critical point at lattice size L can be defined according to finite-size scaling 
theory (Barber 1983) as the coupling xl such that 

Rl{xl) = 1 (3.1) 

where Rl{x) is the scaled energy-gap ratio 

and Fl{x) is the energy gap for lattice size L. This point is found by calculating the energy 
eigenvalues at a cluster of 5 equally spaced points in the neighbourhood of xl, and then 
finding xl by interpolation between them. The spacing between the points was chosen as 
Ax = 0.001, estimated to balance the truncation and round-off errors in the calculation. 
The values of all other observables can then be estimated at x^ by the same finite-difference 
interpolation procedures. Tables 1 and 2 list the pseudo-critical points xl, and the values 
of the calculated observables at coupling xl for each pair of lattice sizes L and (L-1), for L 
= 2,3,4,5 and 6. The values of x^ for L = 2 to 5 listed in Table 1 agree through six figures 
with those calculated previously (Hamer 1983). 

Table 1 lists values for the ground-state energy per site eo,L for lattice size L, and its 
derivatives e'g^ and e'^ where the prime denotes differentiation with respect to x. The 
values are expected to be accurate to the figures quoted (or better) as regards round-off 
error. The truncation error in the 5-point interpolation process is harder to estimate, since 
it involves unknown higher derivatives of eo, but we estimate it should be no more than 
about 1 part in 10^^ for eo and 1 part in 10^ for eg. 

We have also listed values in Table 1 for the magnetic susceptibility, defined by 

Xl - -j2 01^2 l'^=o y'^-'^) 

This derivative was also estimated by a finite difference method, using a cluster of 5 data 
points around h = 0, with a spacing Ah = 0.0003, giving an estimated truncation error of 



4 



no more than 1 part in 10^ in the susceptibihty. This calculation was a little too large to 
carry through for L = 6, with the facilities available. 

Table 2 lists the energy gap between the odd and even sectors, and its derivatives 
and F'l, at each xl- Values are also listed here for the quantity defined by 

Mz. = (0|ai(l)|l) (3.4) 

where |0), |1) are the lowest-lying energy eigenvectors in the even and odd sectors, respec- 
tively. It can be shown (Yang 1952, Uzelac 1980, Hamer 1982) that this quantity converges 
towards the spontaneous magnetization in the bulk limit. Unfortunately, for technical rea- 
sons we were again unable to calculate this quantity for L=6. Since the accuracy of the 
wavefunction is only the square root of that of the eigenvalue, the round-off error in these 
values is expected to be about 1 part in 10^. 



B. Critical Point 



The sequence of pseudo-critical points xl converges rapidly, as can be seen in Fig. 1, 
where xl is plotted against 1/L^. To estimate the bulk limit (L — oo) of this sequence, 
we have employed various algorithms discussed above, as well as a simple polynomial fit in 
and higher powers. 

Our final estimate of the critical point is 

= 0.32841(2) (3.5) 

This is consistent with our earlier finite-size estimate of Xc = 0.3289(10) (Hamer 1983), 
but nearly two orders of magnitude more accurate. Henkel (1987) obtained an improved 
estimate Xc = 0.3282(1) from lattices up to 5x5 sites. 



C. Critical Indices 

Finite-size scaling theory (Barber 1983) also tells us how to estimate the critical indices 
for the model. The finite-lattice susceptibility xl, for instance, is predicted to scale at the 
critical point like 

Xl{xc) ~ L^oo (3.6) 

and hence one finds that 

i^(l - ^^^) ~ --, L^oo (3.7) 
Xl-i[Xl) V 



Similarly, ratios of the finite-lattice 'magnetizations' (equation |3.4| ) give estimates of /3/z/. 
Finally, estimates of the index \ jv can be obtained from the Callan-Symanzik 'beta function' 
(Barber 1983), 
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via 



i-oo (3.9) 

One would expect to obtain estimates of the ratio ajv in a similar fashion from the 
'specific heat', 

^^^^^ = -Z^S^' ^^-'^^ 

but it is known (Hamer 1983) that these estimates are very poor, too high by a factor of 
nearly 2. The reason is easily found: the ground-state energy or specific heat contains a 
'regular' or analytic piece as well as the singular term (Privman and Fisher 1984). Henkel 
(1987) has cleverly sidestepped this problem, using a transition amplitude to find ajv^ in 
analogy to equation (|3.4|) . Here, we eliminate the regular term by subtracting: 



4l - ~ ^"/'^-S ^-oo, (3.11) 

and using successive ratios of these differences to estimate X — ajv The estimates so obtained 
for the critical index ratios are listed in Table 3. 

Alternatively, 'logarithmic' estimates of the critical indices may be obtained as follows: 

hl\XL{xL)IXL-\{xL)\ 7 . _ /oioN 

ln[L/(L-l)] V ^ ' 

These alternative estimates are listed in Table 4. The finite-size corrections are generally 
smaller for the logarithmic estimates. 

These finite-size estimates of the critical indices agree closely with the previous calcula- 
tion of Hamer (1983) up to L = 5; and remarkably enough, most of them agree to within 4 
significant figures with the equivalent results obtained for the triangular lattice (Hamer and 
Johnson 1986, Price et al 1993). 

The same algorithms mentioned above have been employed to extrapolate these se- 
quences to their bulk limit. The sequences are very short, and may have slight irregularities, 
so that the tabular algorithms are generally no more accurate than simple graphical meth- 
ods or polynomial fits in the extrapolation. The resulting estimates are listed at the foot 
of Tables 3 and 4. The errors in these estimates are inevitably rather subjective, but the 
variation between different algorithms gives some indication of the likely error. 

Figure 2 graphs the estimates of 1/z/ from Tables 3 and 4 as a function of 1/L: it 
can be seen that the behaviour is almost precisely linear for the estimates from Table 3. 
Correspondingly, the Neville tables and polynomial fits give stable results, while the Lubkin 
and Bulirsch-Stoer algorithms give less stable results, possibly a little higher. We conclude 
that 

- = 1.591(2) (3.13) 

V 

The estimates for ajv are not quite so well-behaved, but our final estimate is 
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- = 0.16(1) (3.14) 

V 

This is a much better resuh than can be obtained directly from the specific heat, equation 

The estimates for the other indices /?/z/ and ^jv are rapidly convergent, but we only 
have data up to L = 5, which were known previously. We find 

- = 0.522(2) (3.15) 

V 

^ = 1.96(1) (3.16) 



D. Energy Amplitudes 

The finite-size behaviour of the energy gap at the critical point is 

Fl[xl)^^, L-^^ (3.17) 
so the amplitude Ax can be estimated by 

LFi^ixj^) ~ Ai, L ^ oo (3.18) 

The sequence of estimates for Ax is shown in Figure 3. It extrapolates to a value 

Ax = 1.39(1) (3.19) 

A value of 1.42 was previously estimated by Henkel (1987). 

The finite-size scaling behaviour of the ground-state energy per site eo at the pseudo- 
critical point is shown in Figure 4. The finite-size scaling corrections appear to decrease like 
in accordance with the Privman-Fisher scaling hypothesis (1984), which states that 
the singular part of the free energy density of a system of finite size L should scale as L~'^ 
(here d = 3). A polynomial fit on this assumption gives 

eoA^i) ~ ^0 - ^) L ^ oo (3.20) 

with 

eo = -0.624(1) (3.21) 

and 

Aq = 0.38(5) (3.22) 
Further evidence for this power-law behaviour can be obtained as follows. Suppose 

eo,L(a:L) ~e*-Ao/LP, L ^ oo (3.23) 
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then 



L1-- —\=pLr^p as L ^ oo (3.24) 

(eo,L-i(a;L-i) - £o,l-2{xl-i)) 

and 

ln[ ^0-^^^^) - eo,L-iixL) - 1] ^ -p asL^oo (3.25) 

eo,L-l{XL-l) - eo,L_2(XL-i) 

The sequences of finite lattice estimates for p are shown in Figure 5. It can be seen that 
the 'hnear' sequence comes down towards 3 from above, whilst the 'logarithmic' sequence 
comes up towards 3 from below. The sequences are a little irregular, however, and the best 
estimate we can obtain for the bulk limit is 

p = 2.8(2) (3.26) 

a little lower than, but still consistent with 3. 

Assuming that p = 3, the scaling amplitude Aq for the ground-state energy can be found 

by 

— {(^o,l{^l)-^o,l-i{xl))^Ao, L ^ oo (3.27) 

The sequence of estimates for Aq is graphed in Figure 6, and extrapolates to a value 

Ao = 0.35(2), (3.28) 

which is in reasonable agreement with equation ( p.22| ). Henkel (1987) previously obtained 
an estimate of 0.39 for this quantity (allowing for the different normalization of his Hamil- 
tonian). Mon (1985) obtained a Monte Carlo estimate of the corresponding free energy 
amplitude in the 3D classical model. 

In order to 'calibrate' this result, we need to know the "speed of light' v, or in other 
words the scale factor needed in this model to make the long-range correlations isotropic in 
space and time at the critical point. We have attempted to estimate this using the dispersion 
relation for the lowest excited state at the critical point, expected to be of the form 

E{k) = vk (3.29) 

in the bulk system. We have calculated the finite-lattice eigenvalues for low-lying excited 
states with non-zero momntum for lattice sizes L=2 to 5, and set 

L 2tt 

VL = —(^L{xL,—)-FL{xL,0))r^V, L ^ OO (3.30) 

Ztt L 

where Fl{x, k) is the energy at coupling x for momentum k. Figure 7 shows the sequence 
of finite-lattice estimates for v as a function of 1/L. They extrapolate to a bulk value 

V = 0.99(3) (3.31) 
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The ratio Aq/v should be a universal number, independent of the normalization of the 
Hamiltonian. From the results above, we find 

— = 0.35(2), (3.32) 

V 

to be compared with values of 0.719 expected according to effective field theory for a single 
free boson (Hascnfratz and Niedermayer 1993), or 0.211 for a single free fermion degree 
of freedom (Appendix). The result (3.31) matches neither of these values. This is not 
surprising, since the effective field theory at the critical point is expected to be a non-trivial 
interacting theory. It might be possible to estimate this quantity via the e-expansion, using 
a Landau-Ginzburg effective field theory. This has not yet been done, as far as we are aware. 



IV. CONCLUSIONS 

We have calculated the lowest-lying energy eigenvalues of the transverse Ising model on 
the square lattice with periodic boundary conditions for lattice sizes up to 6x6 sites, using 
the conjugate gradient method. Finite-size scaling theory has been employed to estimate 
the critical parameters, which are compared with previous estimates in Table 5. 

It can be seen that our present estimates agree well with earlier finite-size scaling results. 
We have achieved a substantial increase in accuracy for the critical point, but only a more 
modest increase for the critical index u. The results appear very compatible with previous 
series analyses, and also with recent estimates for the classical 3D Ising model, and field 
theory. This provides further confirmation of the universality between these transitions. 
Finally, it can be seen that the accuracy of the exponents for the quantum model is now not 
very far behind that for the classical model. 

We have also estimated the finite-size scaling amplitudes for the energy eigenvalues at 
the critical point. For the spin gap we find 

Ai = 1.39(1) (4.1) 

to be compared with a previous estimate by Henkel (1987) of Ai — 1.42. 
For the ground-state energy, we have shown evidence that 

~ Co - L ^ oo (4.2) 

and estimated 

Ao = 0.35(2), (4.3) 

to be compared with a previous estimate of 0.39 by Henkel (1987). It should be possible to 
predict this amplitude from Landau-Ginzburg effective field theory. 

An extension to 7x7 sites of these exact diagonalization calculations is hardly feasible at 
the present time, but there are some very precise approximate methods now available, such 
as the density matrix renormahzation group (White 1992) and path integral Monte Garlo 
techniques (Sandvik 1992). These might well be able to extend the results to larger lattice 
sizes, and allow much improved finite-size scaling estimates of the critical parameters. They 
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could also confirm whether or not the Casimir energy scales as in equation ( [4.2|) . Our exact 
diagonalization results should provide a useful calibration for such studies. We look forward 
to seeing such calculations in the future. 

We have chosen here to work on the square lattice rather than the triangular one, because 
the Hamiltonian matrix is somewhat smaller, and the leading finite-size corrections are 
expected to be much the same for both lattices. There are, however, some hints of irregularity 
or alternating behaviour in some of the square lattice sequences. It might well be that the 
triangular lattice results are smoother. 
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Appendix 

The finite-size scaling amplitude for the ground-state energy ( "Casimir amplitude" ) can 
be calculated for free fields as follows. 

1. Free boson case 

The zero-point energy of a free boson field is given in d space dimensions by 

^o = ^E^fc (4) 
^ k 

i.e. uJk/'2 for each momentum mode. On a lattice, the free particle Hamiltonian can be 
written in a finite-difference form 

where the lattice spacing has been set to 1. The eigenmodes are plane waves 

0(n, t) = ^ El^'ce^^'^--^*) + ale-^C^-"--^*)] (6) 

k 

whence 
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a;, = [2^(1- COS fc,)]'/' (7) 



i=l 



and for periodic boundary conditions the allowed momenta are 



27r 

ki^—h, 0,1,2,- ■■ (8) 



for a lattice oi N = L''' sites. Hence 

d 



^o = ^E[2E(l-cosfcO]'/' (9) 

k 1=1 



I.e. 



^ {li}=Q i=l ^ 

Now the leading finite-size correction to this sum arises from the infrared (small momentum) 

behaviour of the lattice sum (Hascnfratz and Leutwyler 1990), and does not depend on the 
cutoff or regularization at large momentum. Thus we may approximate for our purposes 

-\ CO d 1 

^o^iz E iUj^n^' (11) 

{li}=-oo 1=1 

Now use the Poisson resummation formula 

+00 1 +OC 2™ 

E n^L) = E 9C-^) (12) 

^ n=-oo ^ 



with 



to show 



gik) = I e^^-V(x)(i'^x (13) 



en = — > / kdk- — '- — rjT^TT — L\m\) (14) 



T^( d+l \ +00 1 

The dash here implies removal of the term m = 0, which corresponds to the (infinite, 
non-universal) bulk ground-state energy per site, which we simply drop. 

The sum involved here is a generalization of the Riemann zeta function. It gives 

4 = -7^ (16) 
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where ior d — 1, Ao is easily evaluated 

Ao^-^ 0.5236 (17) 
6 

the result familiar from conformal field theory. For higher dimensions, we have evaluated 
the sum numerically 

d = 2 : = 0.7189 (18) 

d = 3 : ^0 = 0.8375 (19) 
The result for d = 2 was given previously by Hasenfratz and Niedermayer (1993) . 

2. Free fermion case 

A similar naive argument can be given for the case of a single species of free Weyl 
(spinless) fermions. The filled Dirac sea has energy 

k 

where again 

a;, = [2^(1- cos A;,)]^/^ (21) 

i=l 

and we assume antiperiodic boundary conditions for the fermions 

k, = ^{2k + l), /, = 0,1,2,- •• (22) 

Hence 

^0^-^ E (23) 

{ti}=— OO 1=1 

Use the Poisson resummation formula again to find 

e'o--^ E (-1)^.-' r k'dk (x = L|m|) (24) 



{mi}'=-oo 



For d—lAois easily evaluated to give 



r(^) ^ (-i)^''"' .... 



OO 



TT 



^0 = = 0.2618 (26) 
also familiar from conformal field theory; while for higher dimensions, we find numerically 

d = 2 : ^ = 0.2106 (27) 

d = 3 : v4o = 0.1957 (28) 
These numbers have not been obtained previously, as far as we are aware. 
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FIGURES 

Figure 1. Graph of the finite- lattice pseudo-critical points as a function of 
The fine is merely to guide the eye. 

Figure 2. Finite-lattice estimates of the index 1/z/ graphed against 1/ L. The full circles 
are 'linear' estimates, the open circles are 'logarithmic' estimates. The lines are merely to 
guide the eye. 

Figure 3. The energy gap amplitude A^^l graphed against 

Figure 4. The ground-state energy per site at lattice size L graphed against 

Figure 5. The effective exponent pl graphed against 1/L. 

Figure 6. The Gasimir amplitude Aq^l graphed against 1/L. 

Figure 7. Finite-lattice estimates of the 'speed of light' v graphed against 1/L. 
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TABLES 



TABLE I. Finite-lattice data at the pseudo-critical points xl, calculated for the pair of lattice 
sites L and (L-1) in each case. Given are the ground-state energy per site eo, its first two derivatives 
eg and eg with respect to x, and the susceptibility %. 



XL 


L 


eo 


e'o 


^0 


X 


0.26034238222 


1 


-0.520684764436 


-2.00000000 


0.00000 


-0.09000 




2 


-0.074060535180 


-0.61163273 


-2.90881 


-0.32656 


0.31600008772 


2 


-0.112709188696 


-0.77866508 


-3.06919 


-0.453734 




3 


-0.070818161984 


-0.59879739 


-4.15119 


-1.00929 


0.32424925229 


3 


-0.075901188836 


-0.63384113 


-4.34370 


-1.13177 




4 


-0.066430308096 


-0.57114518 


-5.09931 


-1.99844 


0.32669593806 


4 


-0.067843120636 


-0.58378935 


-5.23641 


-2.11883 




5 


-0.064637823298 


-0.55295527 


-5.83324 


-3.29209 


0.32758326752 


5 


-0.065130784765 


-0.55817035 


-5.92135 






6 


-0.063752757694 


-0.54012509 


-6.42027 





TABLE IL Finite- lattice data as in Table 1, consisting of the mass gap F, its first two 
derivatives^' andF" with respect to x, and the 'magnetization' M. 



XL 


L 


F 


F> 




M 


0.26034238222 


1 


2.000000000000 


0.0000000 


0.00000 


1.00000 




2 


1.000000000000 


-3.4007921 


6.05294 


0.67337 


0.31600008772 


2 


0.820891162135 


-3.0223297 


7.44433 


0.71908 




3 


0.547260774756 


-4.4695972 


13.9564 


0.58357 


0.32424925229 


3 


0.510885529302 


-4.3470817 


15.7416 


0.59705 




4 


0.383164146983 


-5.4088281 


25.0603 


0.51520 


0.32669593806 


4 


0.370007329425 


-5.3452110 


26.9450 


0.52130 




5 


0.296005863592 


-6.2337554 


39.3896 


0.46488 


0.32758326752 


5 


0.290490201648 


-6.1980387 


41.1167 






6 


0.2420751678 


-6.9850090 


56.7636 
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TABLE III. Finite-size scaling estimates of the critical indices, as defined by equation (3.7) 
and following, where L is the larger of the two lattice sizes used in the estimate. 



T -1 / -t 1 

L 1/u 1 — a/u 




7/1/ 


P 


2 1.27817 


0.653264 


1.44880 




3 1.38021 1.88408 


0.565352 


1.65132 


1.71862 


4 1.43242 1.20661 


0.548326 


1.73469 


2.09567 


5 1.46377 1.05069 


0.541083 


1.78195 


2.30781 








2.42047 


oo 1.591(2) 0.84(1) 


0.523(2) 


1.95(1) 


2.8(2) 


lABLb IV. Logarithmic hmte-size scaling 


estimates of the critical indices, as 


defined by 


equation (|3.12|) and following, where L is the lar^ 


5er of the two lattice sizes used in the estimate. 


L 1/z^ 1 — a/z/ 


(i/u 




P 


2 

3 1.52002 2.43901 


0.514989 


1.97178 


4.83686709 


4 1.54105 1.24804 


0.512493 


1.97642 


4.16837815 


5 1.55226 1.05715 


0.513266 


1.97479 


3.85524695 


6 1.55937 0.98287 






3.63001774 


cx) 1.593(3) 0.84(1) 


0.521(3) 


1.96(1) 


2.5(2) 


TABLE V. A comparison of critical parameters obtained in the present work with some 


others obtained elsewhere. Key: HT = high-temperature series; LT = low-temperature series; FS 


= finite-size scaling; MC = Monte Carlo; TR = triangular lattice; SQ = square lattice 


u 


a 






(2+1)- dimensional Ising model 








FS SQ Hamer (1983) 0.635(5) 






0.3289(10) 


FS SQ Henkel (1984, 1987) 0.629(2) 


0.11(1) 


0.324(9) 


0.3282(1) 


FS SQ This work 0.629(1) 


0.10(1) 


0.328(2) 1.23(1) 


0.32841(2) 


FS TR Hamer & Johnson (1986) 0.627(4) 




0.332(6) 1.236(8) 




FS TR Price et al (1993) 0.627(2) 


0.12(2) 


0.324(3) 1.23(1) 




HT SQ He et al (1990) 0.637(4) 


0.11(2) 


1.244(4) 


0.32851(8) 


LT SQ Oitmaa et al (1991) 0.64(3) 


0.096(6) 


0.318(4) 1.25(2) 




3-dimensional Ising model 








HT BCC Butera and Comi (1997) 0.6308(5) 




1.2384(6) 




MC Hasenbusch (1999) 0.6296(3) (4) 




1.2367(11) 




Field Theory 








Guida and Zinn-Justin (1997) 0.6304(13) 


0.109(4) 


0.3258(14) 1.2396(13) 
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